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A concept of space radiation-voltaic batteries (SRVBs) based on various space radiation environments, in- 
cluding the earth radiation belt (ERB), solar proton event (SPE), and galactic cosmic rays (GCRs) is presented 
in this paper. The energy deposition, structural parameter optimization, and output performance of SRVBs are 
studied through the Monte Carlo toolkit Geant4 and theoretical formulas. Energy deposition in the three radia- 
tion environments monotonically decreases as the depth in both silicon film models and energy conversion units 
of the SRVBs increase. The recommended junction depths of SRVBs based on ERB, SPE, and GCR are 0.5, 
3.0, and 0.5 um, respectively. In the proposed design scheme, the maximum output power density of SRVBs in 
ERB, SPE, and GCR are 2.022107’, 1.113x1073, and 3.844107!” W-cm7?, respectively. Calculation results 
suggest that SRVBs based on ERB and SPE are high-potential candidates for space power sources. The results 
could help raise awareness on the use of space radiation particle energy and facilitate research on SRVBs. 
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I. INTRODUCTION 


An electric power source can be created by exposing a 
p-n junction to radioactive particles, such as alpha and beta 
sources. Becker and Kruppke (1937) observed that electron- 
hole pairs (EHPs) are produced in selenium semiconduc- 
tor materials through cathode electron bombardment [1]. 
This phenomenon is called the electron-voltaic effect, which 
Ehrenberg and Billington studied. Rappaport (1953) found 
that isotope decay beta particles also produce EHPs in p-n 
junction devices and called this phenomenon the betavoltaic 
effect [2]. The concept of an alphavoltaic battery was pro- 
posed by Pfann and Roosbroeck in 1954 [3]. Thereafter, 
many studies investigated betavoltaic [4-14] and alphavoltaic 
[15] batteries. 

Radiation sources used in previous studies on nuclear bat- 
teries include cathode electron tubes, artificial or natural ra- 
dioisotopes, and particle accelerators. Research on nuclear 
batteries has mostly focused on solid, liquid, or gaseous ra- 
dioisotopes from the earth. An interbedded epitaxial silicon 
betavoltaic based on an artificial radioisotope Ni was dis- 
cussed through experiments and a MCNP5S simulation [11]. 
The theoretical maximum efficiencies of betavoltaics that uti- 
lize artificial radioisotopes 35S, 9S, and VY were studied 
by employing Monte Carlo codes, including Geant4, Pene- 
lope, and MCNPX [13]. A customized low-energy electron 
accelerator was used to characterize the performance of be- 
tavoltaics by varying the electron beam energy and electron 
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beam flux [14]. In space, radioactive particles, such as pro- 
tons, electrons, and helium particles, emit a large amount of 
energy. These particles have great potential as energy sources 
for nuclear batteries. However, little is known about the use 
of the energy of these space radiation particles. 

This paper presents a concept of space radiation-voltaic 
batteries (SRVBs) based on various space radiation environ- 
ments, including the earth radiation belt (ERB), solar pro- 
ton event (SPE), and galactic cosmic rays (GCRs). SRVBs 
convert the kinetic energy of particles in space directly into 
electric power. This study analyzes the energy deposition, 
structural optimization design, and output performance of 
SRVBs through the Monte Carlo toolkit Geant4 and theoret- 
ical formulas. The results demonstrate that SRVBs based on 
ERB and SPE are high-potential candidates for space power 
sources. 


I. METHODS 


A silicon p-n junction, which has low cost and mature pro- 
cessing technology, was used as the energy conversion unit 
(ECU) of SRVBs. The AE-8/AP-8 model in ERB and the 
data of SPE in 1972 were employed in this study. Due to 
the particles’ high-flux proportion, proton and helium parti- 
cles were simulated in the radiation environment of GCR. The 
GCR model is CREME 86. 

To ensure that the ECU is located in the area of maximum 
energy deposition, silicon film computation models combined 
with various space radiation environments were constructed 
(Fig. la). The transport of proton, electron, and helium par- 
ticles in silicon semiconductor materials was then simulated 
with the use of the Monte Carlo toolkit Geant4 in the three 
space radiation environments. Based on the simulation re- 
sults, SRVB computation models (Fig. 1b) were constructed. 


S010402-1 


LIU Yun-Peng et al. 


7\\ ry 
Space Radiation 
Environments ral 


Energy 
deposition 
layer with 
0.02 cm thick 


Circular Silicon Film 


Cross area is 1 cm? with 3 em thick 


(a) 


FIN ry 


Space Radiation 
Environments 


Energy 

deposition 

layer with Circular Si-based ECU 

0.5 pm thick : ra f 
Cross area is 1 cm? with 1000 pm thick 


(b) 


Fig. 1. (Color online) (a) Silicon film and (b) SRVB computation 
models in ERB, SPE, and GCR. 
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where q indicates the quantity of electric charge of an elec- 
tron or a hole (in C), S indicates the cross section of the ECU 
(in cm?), Eehp is the average energy required to generate an 
EHP (in MeV; the Eenp of Si is 3.640x 10-3 MeV), E(n) is the 
energy deposited onto the n” layer of the silicon ECU by par- 
ticles in the three space radiation environments (in MeV), the 
subscript k is the number of types of particles, and @,; is the 
total flux of the k” type of particles (in cm~*s7!), In ERB, the 
fluxes of electron and proton particles are 4.648010’ and 
6.3873x10’ cm~s"!, respectively. In SPE, the flux of pro- 
ton particles is 2.4506x10!° cm~*s~!. In GCR, the fluxes of 
proton and helium particles are 2.916 and 1.131 cm~?s7!, re- 
spectively. 


With the use of MATLAB, other parameters of the SRVB 
were obtained according to empirical formulas. These pa- 
rameters include the optimized doping concentration of the p- 
type (N4) and n-type regions (Np), open-circuit voltage (Voc), 
maximum power density (Pn), fill factor (FF), and conver- 
sion efficiency (77). The calculations made in the current study 
are similar to those in our other papers [9-11]. 


Nucl. Sci. Tech. 25, S010402 (2014) 


The study assumed that all energy deposited onto the ECU 
was used to produce EHPs despite the ECU radiation damage 
caused by high-energy particles. 

In SRVB computation models, the 1 000 um thick Si-based 
ECU was segmented into 2 000 energy deposition layers of 
thickness 0.5 um. One or several energy deposition layers in 
the incident surface of the Si-based ECU represent a heavily 
doped region, that is, the junction depth. Moreover, the other 
energy deposition layers represent a lightly doped substrate 
region. 

The radiation-generated current density of the SRVB, that 
is, its short-circuit current density (Jsc), was calculated by 
collecting EHPs from the energy deposition distribution on 
each layer of the ECUs based on ERB, SPE, and GCR. 
The collection ratio equation of EHPs can be written as fol- 
lows [16]: 


CE (n) = 1 — tanh (x,/L), (1) 


where CE(n) is the collection ratio of the EHPs of the nth 
energy deposition layer in the SRVB computation model, L 
is the minority carrier diffusion length (in um), and x, is the 
distance between the position of the nth layer of Si and the 
depletion region (in um). If the collection ratio of the elec- 
tron hole pairs in the depletion region is 100%, Jsc can be 
expressed as follows: 


X, CE (n) x (E1 (n) x By + Ex (n) X By + +++ + Ex (n) x Bp), (2) 


Ill. RESULTS AND DISCUSSION 


Figure 2 shows that the silicon film irradiated by ERB is the 
fastest one to achieve maximum accumulated energy deposi- 
tion, whereas the film irradiated by GCR is the slowest one. 
Many particles in ERB are distributed primarily in the low- 
energy section, while many particles in GCR are distributed 
principally in the high-energy section. The energy spectrum 
of ERB is narrower than that of GCR. Thus, energy from ERB 
is more easily and quickly deposited onto the shallow part of 
the silicon film than that from GCR. However, energy deposi- 
tion monotonically decreases as the depth in the silicon film in 
all space radiation environments increases. Consequently, the 
ECUs were designed to accept irradiation directly from ERB, 
SPE, and GCR (Fig. 1b). Energy deposition distribution in 
the ECU based on ERB, SPE, and GCR is obtained accord- 
ing to the SRVB computation model (Fig. 3). The trend of 
energy deposition distribution in the ECU is similar to that in 
the 3 cm-thick silicon film, as expected. 


The conversion efficiency of the nuclear battery is associ- 
ated with the flux of radiation particles. Generally, a high flux 
of particles increases battery efficiency. Optimization results 
(Table 1) show that the SRVBs in ERB and SPE have high 
maximum output power density. However, the maximum 


S010402-2 


ENERGY DEPOSITION, PARAMETER... 


Nucl. Sci. Tech. 25, S010402 (2014) 


Table 1. Optimization output performance and corresponding parameter values of SRVBs irradiated by ERB, SPE, and GCR at different 


junction depths at 300 K 
eee Saas Nacom®) ——Np(em3)—Jgc(Acm™2)—Voc(V)  PaWem?)  FF(%) n% 
0.5 4.467x10!8 1.000x10!6 7.299x1077 0.3664 2.022 1077 75.61 6.850 
1.0 1.955x10!8 8.913x10!5 7.320 1077 0.3642 2.013107? 75.51 6.819 
ERB 1.5 1.259x10!8 8.913x10!5 7.311x1077 0.3634 2.005x1077 75.47 6.792 
2.0 8.913x10!” 7.943x10!5 7.333x1077 0.3615 1.998x1077 75.38 6.768 
2.5 6.31010!” 7.943x10!5 7.333x1077 0.3605 1.992x1077 75.34 6.747 
3.0 5.623x10!7 7.943x10!5 7322x1077 0.3602 1.986x1077 75.32 6.727 
0.5 1.259x10!9 8.913x10!3 2.649x 1073 0.5191 1.112x1073 80.86 8.947 
1.0 5.623x10!8 8.913x10!3 2.650x1073 0.5190 1.112x1073 80.86 8.949 
SPE 1.5 3.548x10!8 8.913x10!3 2.651x1073 0.5190 1.112x1073 80.86 8.950 
2.0 2.512x«10!8 8.913x10!3 2.651x107 0.5190 1.113x1073 80.86 8.952 
25 1.778x10!8 8.913x10!3 2.652x107 0.5190 1.113x107 80.85 8.954 
3.0 1.413x10!8 8.913x10!3 2.653x1073 0.5189 1.113x1073 80.85 8.956 
0.5 1.000x 102 1.122x10!5 8.659x107!4 0.001523 3.844 10717 29.16 6.713x1073 
1.0 7.079x10!9 1.122x10!5 8.660x107!4 0.001522 3.843107!” 29.16 6.711x1073 
GCR 1.5 4.467x10!° 1.122x10!5 8.662x107!4 0.001521 3.842x 107!” 29.16 6.709x1073 
2.0 3.162x10!9 1.122x10!5 8.665x107!4 0.001520 3.841x107!7 29.17 6.707x1073 
2.5 2.512x10!° 1.122x10!5 8.667x107!4 0.001519 3.840x10717 29.17 6.706x107 
3.0 1.995x10!° 1.122x10!5 8.670x107!4 0.001518 3.839x107!7 29.17 6.704 1073 
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Fig. 2. (Color online) Energy deposition behavior in silicon films of 
(a) ERB, (b) SPE, and (c) GCR. 
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Fig. 3. (Color online) Energy deposition distribution in ECUs of 
ERB, SPE, and GCR. 


power density and conversion efficiency of the SRVB based 
on GCR are only 3.844x107!7 W/cm? and 6.713x107°%, re- 
spectively. The results show that the energy of GCR is dif- 
ficult to collect and that, by contrast, SRVBs based on ERB 
and SPE are high-potential space power sources. In addition, 
the efficiencies of SRVBs based on ERB and GCR slowly de- 
crease as the junction depth increases, whereas that of SRVB 
based on SPE gradually increases. This increase is associated 
with the energy deposition distribution in the ECU. This result 
suggests that the junction depth of SRVBs in ERB and GCR 
should be shallow and that of SRVB irradiated by SPE should 
be deep. The recommended junction depths of SRVBs based 
on ERB, SPE, and GCR are 0.5, 3.0, and 0.5 um, respectively. 


In space, radiation particles may have low or high ener- 
gies. These particles run in all directions. Therefore, thin 
ECUs have difficulties in absorbing all the energy of radia- 
tion particles. Table 2 shows that only a small amount of en- 
ergy emitted by the space radiation environments, especially 
by GCR, is deposited onto the ECUs. Energy is mostly es- 
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Table 2. Power distribution in three space radiation environments 


Space radiation environments ERB 

Total power (W) 1.222x107> 
Deposition in ECU (W) 2.952x10~° 
Ratio to deposition in ECU (%) 24.17 


caped from the other angle of the ECU, penetrated the ECU, 
or scattered from the surface of the ECU. Adding the num- 
ber of ECUs and optimizing the structure of the ECU sur- 
face ensures that all energy from the space radiation environ- 
ment could be deposited onto the ECU. Under these condi- 
tions, the maximum power density of SRVBs with the rec- 
ommended junction depths based on ERB, SPE, and GCR are 
8.371x1077, 1.865x107?, and 4.437x107!! W/cm?, respec- 
tively. 


IV. CONCLUSIONS 


In this paper, we propose SRVBs based on ERB, SPE, and 
GCR. The optimization output performance and correspond- 
ing parameter values of SRVBs at different junction depths 
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